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Summary 

This report describes a formal method to quantify structural 
damage tolerance and reliability in the presence of a multitude 
of uncertainties in turbine engine components. The method is 
based at the material behavior level where primitive variables 
with their respective scatter ranges are used to describe behav- 
ior. Computational simulation is then used to propagate the 
uncertainties to the structural scale where damage tolerance 
and reliability are usually specified. Several sample cases are 
described to illustrate the effectiveness, versatility, and matu- 
rity of the method. Typical results from this method demon- 
strate that it is mature and that it can be used to probabilistically 
evaluate turbine engine structural components. It may be 
inferred from the results that the method is suitable for 
probabilistically predicting the remaining life in aging or dete- 
riorating structures, for making strategic projections and plans, 
and for achieving better, cheaper, faster products that give 
competitive advantages in world markets. 

Introduction 

Achieving and retaining competitive advantages in world 
markets necessarily leads to proactive drives for getting better, 
cheaper, faster products to market, a phenomenon that becomes 
even more important in the high-tech sector including aero- 
space vehicles. The awareness for natural resource conserva- 
tion also leads to the cost-effective useful life extension of 
existing products. These activities require that we effectively 
use available resources and that we formulate methods to 
quantify the current strength of a specific structure or compo- 
nent and reliably evaluate its remaining strength and/or life. A 
multitude of uncertainties must be dealt with to meet these 
requirements: new, unproven methods for design and manufac- 
turing; a lack of sufficient data on new candidate materials; 
unknown assumptions and conditions related to the initial 
design; records of environmental exposure; and material deg- 
radation from various factors. 

With respect to the aforementioned uncertainties, probabi- 
listic methods offer formal approaches to quantify those uncer- 


tainties and to evaluate their effects on material behavior, on 
service, and on attendant reliabilities and risks. The objective of 
this report is to describe one probabilistic method used to 
evaluate the effects of uncertainties on structural damage 
tolerance and reliability from material behavior to service life 
and to retirement for cause. The probabilistic method is based 
on formulations that describe the physics in terms of primitive 
variables and respective scatter ranges at the lowest engineer- 
ing manageable scale and at all subsequent scales up to the 
highest where reliability is evaluated. The method relies on 
computational simulation to propagate uncertainties from the 
elementary scale and all intermediate scales where the proba- 
bilistic evaluation and reliability of specific responses are 
needed. The method has evolved over the past 1 5 years and has 
matured sufficiently to evaluate structural damage tolerance 
and reliability under various scenarios (refs. 1 to 3). It has 
several unique features, two of which are the most useful for 
presenting results: (1) quantifiable reliability in terms of cumu- 
lative distribution functions and (2) sensitivity factors of the 
primitive variables that affect that reliability. 

This report begins with an introduction of the fundamental 
concept and computational simulation methods and gives a 
simple example. Next is a brief description of the method and 
its attendant computer codes. Several sample cases are then 
discussed to (1) illustrate its versatility, (2) present the large 
amounts of information generated, (3) show the relevance of 
the information, and (4) make recommendations concerning 
design, material development quality, certification, in-service 
health monitoring, retirement for cause, and recycling. The 
description is limited to typical results obtained and their 
respective interpretations. References are cited for specific 
details. 


Fundamental Concepts 

The following simple example describes some fundamental 
concepts of probabilistic structural analysis and design: the 
probabilistic evaluation of the tip displacement of a cantilever 
beam loaded at the free end as shown schematically in 
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figure 1(a). The equation (deterministic model) for predicting 
the tip displacement is 

4Pl~ 


where w is the response variable, P is the load, l is the length, 
E is the material stiffness, b is the width, and t is the thickness. 
This equation describes the physics of the response sought and 
includes the fundamental parameters (primitive variables) that 
govern the tip displacement. These primitive variables can also 
be grouped in three generic categories: load (P), geometry (/, b, 
and 0. material (E). If we make several cantilever specimens, 
there will probably be some scatter of values for each of the 
primitive variables. Therefore, the task of probabilistic simula- 
tion is to account for the effects of that scatter on the displace- 
ment of the beam. 

The task is considerably simplified when we recognize that 
(1) the tip displacement equation is the analog of the physical 
testing machine and (2) the scatter range in the primitive 
variables, P, /, b, t, and E, can be assumed to be represented 


by simple, well-known statistical distributions (fig. 1 (b)), which 
will be helpful in simulating the effects of scatter on beam 
displacement. The following procedure is used to evaluate the 
effects of the uncertainties of the primitive variables on the tip 
displacement: 

1 . Decide the range of the scatter in each primitive variable. 
This range in practical cases is established from experience but 
for this example, assume that scatter is about ±5 percent from 
the assumed mean value. The scatter for the modulus is between 
24 and 28 mpsi; for the length, between 19 and 21 in.; for the 
width, between 0.95 and 1.05 in.; for the thickness, 1.20 and 
1 .30 in.; and for the load, between 80 and 1 20 lb. It is important 
to note that the only test data available were the mean values for 
the primitive variables. The range of the scatter was assumed 
for reasons that will become clear later. Note that the mean 
value for each primitive variable with a truncated distribution 
generally is where the vertical line (drawn from the peak of the 
respective distribution) intersects the horizontal line. 

2. Randomly select for each primitive variable in the equa- 
tion a value from within its respective distribution. Having 
the simple statistical distributions allows nonbiased random 
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Fiqure 1 -Probabilistic structural analysis and response, (a) Structural analysis model, (b) Distribution of 
primitive variables P, l, E,t,b. (c) Distribution of response variables w. (d) Sensitivity factors. 
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selections to be made. For example, the values selected 
randomly can be 25.5 mpsi for E\ 20.8 in. for /; 0.99 in. for b\ 
1.27 in. for t\ and 1 15 lb for P, 

3. Substitute the randomly selected values in the equation to 
obtain 0.08 in. for the tip displacement. 

4. Repeat step 3 for different sets of primitive variable values 
until sufficient data (as described subsequently) have been 
accumulated to plot the probability distribution curve 
(fig. 1(c)). For example, the mean value will be close to 

0.065 in. There is about a 95-percent probability from the 
cumulative probability curve (fig. 1(c)) that the tip displace- 
ment will be less than 0.08 in., or 95 of 1 00 trial calculations in 
step 3. 

The generation of the data as described in step 4 is called a 
direct Monte Carlo simulation and generally requires a large 
number of simulations (tens of thousands). Methods and algo- 
rithms have been developed to generate the two probability 
graphs for the displacement with relatively few simulations, 
which for the first-order reliability method (FORM) is 2/1+1, 
where n is the number of primiti ve variables with scatter ranges. 
One such method, the fast probability integration (FPI, ref. 3), 
was used to generate the probability curve of figure 1(c). The 
application of FPI requires inputs of the mean value, scatter 
range, and the known or assumed probability density function 
of the scatter for each participating primitive variable. It 
becomes evident that the probabilistic simulation can be per- 
formed with known mean values, judiciously assumed scatter 
ranges, and respective distributions for the primitive variables. 

A byproduct of the FPI is the sensitivity factors (fig. 1(d)). 
These factors probabilistically quantify and order the sensitiv- 
ity of the cumulative distribution function of the response 
variable (displacement) to the uncertainty (scatter range) in the 
primitive variables (material, geometry, loads) . For this example, 
the load (primitive variable) has about the same effect on the tip 
displacement (response variable) as the geometry parameters 
(primitive variables) at a low probability (<1 in 1000) whereas 
the thickness (primitive variable) dominates at high probabili- 
ties (>999 in 1000). Sensitivities are discussed in the following 
sections. For application to structural components or systems, 
the foregoing method using FPI is generalized as follows: 

1. Develop or use a deterministic model for the entire com- 
ponent or system with its boundary, load, and expected envi- 
ronmental conditions. Practical structural situations would 
dictate that this be predominately a finite-element model. 

2. Identify the primitive variables in the deterministic model. 
These will include material properties, fabrication process 
variables (pressure, temperature, and other loading condi- 
tions), structural parameters, loads (including environment), 
boundary conditions, and so forth. For composite structures, 
use integrated composite mechanics to predict the composite 
properties, beginning with micromechanics and accounting for 
both fabrication and environmental effects. 


3. Obtain or assume mean values, scatter range, and proba- 
bilistic distribution (density function and standard 
deviations) for each primitive variable. 

4. Perturb each primitive variable on either side of its respec- 
tive means by a reasonably small amount (usually up to 
1 0 percent) as a ratio of the respective standard deviation (up to 
20 percent may be used, but with caution). Any amount greater 
than 20 percent may be more a shift or even multimodal instead 
of reasonable scatter. Amounts greater than 20 percent would 
indicate that the input data might represent more than one 
population. 

5. Conduct deterministic analyses with the values selected in 
step 4. 

6. Repeat steps 4 and 5 several times (2n + 1 for FORM) to 
generate sufficient information for FPI use. 

7. Use FPI to generate the probability distribution functions 
for the desired responses, displacement, stress, and frequency 
and for the respective sensitivities at select probability levels. 
Recall that the number of perturbations usually needed with FPI 
is 2n+l, where n is the number of primitive variables and 1 is 
using only their means in the same simulation run. 

The above generalized method is practical through the use of 
computer codes to be described subsequently, and it is appli- 
cable to almost all disciplines and structures described herein. 

Probabilistic Simulation of Components 
and System Reliability 

To evaluate turbine engine component and system durability 
and reliability, probabilistic simulations must be performed for 
the loading conditions, the structural components, including 
supports, and the material(s) behavior. Each of the 
simulations must be defined by inputs of its respective deter- 
ministic model, primitive variables, and their attendant scatter 
range. The probabilistic simulation proceeds to evaluate a 
specific response and its scatter. The evaluated response is then 
compared with the corresponding resistance to assess the 
probability of failure, which can be used later to evaluate 
component and system durability and reliability. Figure 2, a 
conceptual schematic of the procedure, shows the three essen- 
tial parts of component and system reliability simulation; the 
structural response obtained; the resistance evaluated; the 
probable damage (overlap of response scatter with resistance 
scatter); the information passed on for reliability and risk 
assessments; and the institutions that participated to develop 
the requisite formalism and then implemented it in an opera- 
tional computational procedure (ref. 1). A block diagram of 
the computer code logic is shown in figure 3. 

The schematics in figures 2 and 3 summarize probabilistic 
structural performance assessment. The concept is relatively 
straightforward and appears simple; however, implementation 
in a workable computer code requires knowledge of advanced 
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structural mechanics, efficient probabilistic algorithms, mate- 
rial behavior, and proficient and subtle computer programming 
techniques. 

In figure 3 , note that uncertainties in the human factor and the 
computer code can also be included; inputs for required perfor- 
mance, component and system longevity, and acceptable reli- 
ability and risk must be provided; and the simulation provides 
information to probabilistically select verification tests to 
assure component and system certification with an acceptable 
reliability and an affordable risk. These will now be described. 

Simple loading conditions can be input directly to the proba- 
bilistic structural analysis. Complex loading conditions require 
system-specific computer codes. Those for the space shuttle 
main engine are simulated by the composite load spectra 
computer code (ref. 2). Probabilistic structural analysis is 
performed by a specialty computer code (refs. 2 and 3). 

Probabilistic Simulation of Material 
Behavior 

Developed at the Glenn Research Center, the probabilistic 
simulation of material behavior is relatively new and as far as 
the author knows, is the only one of its kind. Since the subject 


of this report is durability and reliability, materials-based life 
prediction is an important part. Therefore, the probabilistic 
material behavior models (PMBM) used in the simulation will 
be described in detail. The deterministic model evolved during 
research on high-temperature metal matrix composites (ref. 4). 
Implementing the deterministic model for probabilistic simula- 
tions was funded from a research grant with the objective 
of formally describing uncertainties in material behavior for 
space shuttle main engine components (ref. 5). 

Conceptually, the model is based on the self-evident axiom 
that each material characteristic property observed by conven- 
tional testing constitutes a multidimensional surface. That 
surface is described by a multidimensional vector for which 
each component represents one observed or assumed effect 
(temperature, stress, time, etc.) on that material characteristic 
property. The surface is represented by a respective multifectpr 
model of product form. The product form is assumed^ to 
conveniently represent mutual interactions among the various 
factors. Each factor consists of four different variables or 
parameters: (1) a terminal or final value, (2) a reference value, 
(3) a current value, and (4) an exponent. The exponent is 
selected to represent continuous monatomic behavior so that 
the factor equals unity when the current value equals the 
reference value and approaches either zero or infinity (depend- 
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Figure 3. — Probabilistic structural simulation for assured certification. 


ing on the specific behavior of the factors) when the current and the subscripts are o , reference; F, final; M, mechanical; 

value approaches the final value. A schematic for a unit cell of and T, thermal. The over dot signifies the rate, 

material is shown in figure 4 for composites and homogeneous 
materials. In equation form, the model is represented as 
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P property 

T temperature 

S strength 

R metallurgical reaction 

N number of cycles 

t time 

a stress 



Figure 4. — Concept of unit cell of material behavior through 
multifactor Interaction model (MFIM). 
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To summarize, the multifactor interaction model provides the 
following: 

• Gradual effects during most of the range, rapidly 
degrading near final stages 

• Representation of the in situ behavior for fiber, matrix, 
interphase, and coating 

• Introduction of primitive variables 

• Consistent representation of all in situ constituent proper- 
ties in terms of primitive variables 

• Room-temperature values for reference properties 

• Continuous interphase growth 

• Simultaneous interaction of all primitive variables 

• Adaptability to new materials 

• Amenability to verification inclusive of all properties 

• Adaptability to incremental computational simulation 

Probabilistic results from the model (ref. 5) are shown in 
figure 5 where the cumulative distribution function curves for 
lifetime strengths are given for three temperatures. The curves 
shift to the left and their scatter range increases slightly with 
increasing temperature, as physically would be expected. One 
can infer from these curves that the MFIM could be used in 
conjunction with selective testing to substantially reduce the 
number of tests and the amount of time required to characterize 
material behavior in complex environments. Note that the 



Lifetime strength/reference strength, S/S 0 

Figure 5.— Probabilistic material behavior model 
(PMBM) simulated lifetime strength for nickel-based 
superalloy subjected to 31 62 stress cycles and 
100 hr of creep. 


MFIM is not restricted to the use just described but is generic 
in that each factor can be further substructured to another set of 
factors which may influence a specific factor (i.e., alloying 
elements, processing conditions, tolerance, assembly misfits, 
etc.). This generic use is discussed in reference 6, which 
describes its application to simulate the effects of the human 
factor in structural analysis. Analogously, the MFIM can be 
used for evaluating aging effects on material deterioration 
simply by including deterioration factors. 


Demonstration Cases 

Two-Stage Rotor 

This case demonstrates one direct way to evaluate compo- 
nent and system reliability under multiple failure modes. A 
schematic of the rotor with a summary of the results is shown 
in figure 6. The details of this c ase are describ ed in reference 7, 
but this report will discuss the signific ance of the results and the 
respective sensitivity factors. Four failure modes were evalu- 
ated as noted in figure 6. The survival probability of the rotor 
for each failure mode and the combination of failure modes are 
determined from a special plot of survival probability versus 
remaining resistance. This plot graphically depicts the critical 
failure mode and the system failure mode. As can be seen, the 
critical failure mode is fracture at the rim in 10 000 cycles, 
which coincides with the system failure mode. Also seen is that 
when the burst failure mode has 100 percent survival probabil- 
ity, the system has only about 65 percent. Table I lists the 
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Figure 6— Structural system reliability of two-stage 
rotor with multiple failure modes, (a) Rotor system, 
(b) Critical failure mode and system failure mode. 
All dimensions are in inches. 
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Means ratio = Remaining resistance/initial resistance 


Failure mode 

Applied stress 

Resistance 

1 . Disk burst 

2. Fracture at bore 

3. Fracture at rim 

4. Progressive 
damage 

Average stress 
Fracture life 
Fracture life 
Yielding of ring a 

Burst strength 
10 000 cycles 
10 000 cycles 
Yield strength 


a Yielding of ring affects all other modes through 
mutual interaction. 


Figure 6. — Concluded, (b) Critical failure mode and 
system failure mode. All dimensions are in inches. 


TABLE I,— ROTOR SYSTEM 
RELIABILITY SENSITIVITY 
FACTORS AT 1/1000 
PROBABILITY 


SPEED 

0.850827 

ROTOR DENS 

.438499 

TEMFE 

.170793 

C 

.133702 

RINGY 

.073086 

RK1C 

.061872 

Kt 

.060917 

AO 

.057976 

E_ROT 

.016011 

BURST 

.011983 

A_LCF 

.005132 

EJUN 

.002698 

Nl 

.000386 

RING DENS 

.000008 

TOLER 

0 


sensitivity factors that influence system failure. The left col- 
umn gives the primitive variables included in the evaluations, 
and the right column gives the magnitude of each factor’s 
contribution to the system probability of failure or reliability. 
These sensitivity magnitudes are part of the probabilistic simu- 
lation from using FPI and are evaluated simultaneously with the 
probability, as mentioned in the cantilever example in the 
section Fundamental Concepts. One may observe from this 


table that the four most dominant factors (in order of magni- 
tude) are rotor speed, density, temperature, and crack growth 
constant C. The remaining factors make negligible contribu- 
tions. For example, critical parameters in damage tolerance 
evaluations (initial crack size A 0 and constants C in the fracture 
mode, N i and K f ) are insignificant in the rotor reliability 
assessment. The only critical material property is the rotor 
density . Therefore, recommendations for rotor material suppli- 
ers are that they control the scatter of the rotor material density 
and the thermal expansion coefficient. Also, rotor designers 
must assure that the rotor does not over- speed, that expected 
temperature spikes be accounted for, and that unexpected hot 
spots be avoided. This probabilistic evaluation of the rotor 
illustrates that probabilistic results can be used to provide 
guidelines for material quality control and design consider- 
ations, both of which are essential in product safety, reliability, 
and cost reduction. 


Combustor Liner 

The engine combustor liner to be described could be a part 
for a supersonic aircraft engine. As a result of the combustion 
process, the liner is subjected primarily to thermal loads and to 
some pressure. Figure 7 shows the finite-element model of the 
liner and the thermal loading profile along the combustor. The 
details of this evaluation are described in reference 8. For this 
discussion, assume that the liner is constructed of a cross-ply 
(woven fabric) ceramic matrix composite and that it is designed 
for the avoidance of vibration frequencies and buckling. The 
probabilistic results are presented and the usage of these results 
for damage tolerance and reliability of the combustor liner 
are described. Guidelines and recommendations for material 
selection may be obtained from the sensitivity factors. 

The cumulative distribution function of the first (lowest) 
vibration frequency is shown in figure 8(a). This frequency 
has a mean of 320 cycles/sec and a scatter range from 290 to 
350 Hz, about 60 Hz. The sensitivity factors for two levels of 
probability are shown in figure 8(b). Evident from figure 8(b) 
is that liner material density and shear modulus have a signifi- 
cant effect on the liner frequency; the liner thickness has the 
dominant effect and is even more dominant at higher probabil- 
ity values; and the order of the sensitivity factors is the same at 
low and high probability values. 

The cumulative distribution function for buckling of the 
liner is shown in figure 9(a) and the respective sensitivities are 
shown in figure 9(b). The mean value of the buckling pressure 
is six times that of the operating pressure, or 60 psi. The 
attendant scatter range is from about 45 to 75 psi. The reliability 
of the liner for buckling is 100 percent with no risk since the 
applied mean pressure is only 10 psi, or 6.5a away from 
probable overlap. The factors that exert the greatest influence 
on the buckling load are the liner thickness (geometry vari- 
able), thermal expansion coefficient (material variable), and 
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1300 °F 
(thermal 
load 1) 


2150 °F 
(thermal 
load 2) 


Figure 7.— Finite-element model for aircraft engine 
combustor liner. 


thermal load (loading conditions of 2150 °F). The material 
moduli (axial, hoop, and shear) have a negligible effect. Note 
that the order of the sensitivities remains the same for low and 
high probability values, which means that the buckling load 
probability is linear throughout the scatter range of each primi- 
tive variable. On the basis of these results, the recommendation 
for material suppliers is that they control the thermal expansion 
coefficient. For designers, it is recommended that they control 
the temperature, select the liner thickness to assure that it will 
survive at least 4.5 times the operating pressure, and specify 
proof test pressures of at least 7.5 times the operating pressure 
to assure that the liner will buckle in the first test to verify the 
simulation results. 



1.00 



0.001 0.999 


Probability level 

Figure 8. — Combustor liner probabilistic vibration fre- 
quency and sensitivities, (a) Cumulative distribution 
function of frequency, (b) Sensitivity factors. 


Space Shuttle Main Engine (SSME) High-Pressure Blades 

The space shuttle main engine blades are in the liquid 
hydrogen pump, are relatively small, rotate at about 40 rpm, and 
are subjected cyclically to very cold and very high temperatures 
(thermomechanical fatigue). The blade airfoil with its respec- 
tive operating loading conditions and finite-element model are 
shown in figure 10. The blade has relatively steep spanwise 
thermal and pressure gradients. 

The cyclic temperature and load effects on the blade materi- 
als were simulated by the MFIM described previously (fig. 4). 
The specific values for the various factors used are listed in 
table II. Note that four factors were sufficient for that simula- 
tion: (1) the temperature dependence with an exponent of 1/2^ 
(2) the stress dependence with exponent rt\ (3) the pressure 
cyclic load with exponent/?; and (4) the thermal cyclic load with 
exponent q. The temperature effects exponent was assumed to 
be a constant based on previous studies whereas exponents n, 
/?, and q were assumed to have the scatter shown in table n. 

The damage propagation path caused by 100 000 fatigue 
cycles is shown in figure 11 for two probability levels 
(1/100 000 and 2/10 000). Obviously, the path with the highest 
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0.001 0.999 


Probability level 

Figure 9. — Combustor liner probabilistic buckling and 
sensitivity, (a) Cumulative distribution function for 
buckling, (b) Sensitivity factors. 



Temperature, °F 




probability will most likely occur first. It is noted that several 
other paths are probable with probability levels between the 
two shown in figure 1 1 . However, none was found with a higher 
probability than (2/10 000). The durability, or damage toler- 
ance, of the blade in its operating environment can be simulated 
by using progressive structural fracture (ref. 9). This approach 
requires neither fine finite-element meshes nor traditional 
fracture toughness parameters. The results for the strain energy 
release rate versus the damage state are plotted in figure 12. The 
two major points worth noting in figure 12 are that the damage 
is stable and progresses rather slowly up to damage state 3, and 
the damage progression increases very rapidly from damage 
state 3 to damage state 4. The plot in this figure displays several 
important aspects of structural durability and/or damage toler- 
ance: 

• The blade is damage tolerant up to damage state 3. 

• With continuing operation, the blade will fracture 
(disintegrate) just past damage state 4. 

• The safe design of the blade with 100-percent reliability is 
up to damage state 2. 



Centrifugal loads 


Figure 10. — Finite-element model of space shuttle 
main engine blade with thermomechanical loads. 
Rotational speed, n, 40 000 rpm. 
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TABLE n— MULTIFACTOR INTERACTION MODEL VALUES USED IN 


PROBABILISTIC SIMULATION OF MATERIAL BEHAVIOR OF HIGH- 

PRESSURE_BLADEOF^PACE^HUrn^E^1AINENGINE 

Variable I Distribution I Mean I Standard deviation 



type 


Value 

Percent 

of 

mean 

Temperature, °F 


2750 

51.4 

2.0 

Final, 7V 

Normal 

Initial, T 0 , 

Normal 

68 

2.04 

3.0 

Final strength, 5^, ksi 

Normal 

212.0 

10.6 

5.0 

Initial stress, a 0 

Constant 

0 

0 

0 

Number of mechanical cycles 
Final, N MF 
Initial, 

Lognormal 

10* 

5x10* 

5.0 

Lognormal 

10 3 

50 

5.0 

Number of primitive vari- 

Normal 

0.25 


3.0 

ables with scatter ranges, n 



Exponent for pressure cyclic 

Normal 

.25 


3.0 

load, p 



Exponent for thermal cyclic 
load, q 

Normal 

.25 


3.0 



Probability that path A occurs, 0.00001 Probability that path B occurs, 0.0002 


Figure 11 . Structural reliability of space shuttle main engine blade for 1 000 000 fatigue cycles and probable 

propagation paths to structural fracture. Blade geometry: 0.75 in. long by 1 .0 in. wide at base. 
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Total cost, units Probability of damage initiation 


State Damage 

0 Undamaged 



Damage state 

Figure 1 2.—Damage tolerance of space shuttle main 
engine blade along most probable progressive 
damage path leading to structural fracture. 



Figure 1 3.— Probabilistic risk-cost assessment of space shuttle main engine blade. 
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Figure 14.— System reliability using fault tree simulation. Failure occurs If 
[g 1 < 0 and g 2 < 0] or [g 3 < 0] or [g 4 < 0]. Finite-element model used for 

all bottom events. 


TABLE m —PRIMITIVE VARIABLES USED IN FAULT TREE SIMULATION 
FOR SYSTEM RELIABILITY 



Mean 1 

Standard 

deviation 

Distribution 

Bottom 

event 

Material orientation 




All 

e. 

0.05236 

0.067544 

Normal 

-.03491 

.067544 

Normal 

All 

e. 

’.08727 

.067544 

Normal 

All 

Modulus 




All 

Elastic, E c z 

18.36x10' 

.44595x1 0 6 

Normal 

Shear, Go 2 

18.63x10' 

.4657x1 0 6 

Normal 

All 

Poisson’s ratio, v 0 8 

.386 

.00965 

Normal 

All 

Creep equation coefficient, B 0 

86.0 

.086 

Normal 

Creep 

Density, p 

.805x10"’ 

.493x10 s 

Normal 

Frequency 


‘Defined at room temperature. 


• The blade should be inspected for damage states 
1 and 2. 

• At damage state 2, the blade must be replaced (retired for 
cause) to ensure safe operation of the SSME. 

• Rather than inspect the blade, a costly and time- 
consuming task, in-service health should be monitored 
based on changes in select blade responses (e.g., vibration 
frequencies and vibration mode shapes) to indicate the 
damage state. 

An important observation from the previous discussion is 
that an abundance of information generated by probabilistic 
computational simulation can be judiciously used from con- 
ceptual design to retirement for cause (from cradle to grave). 
Also, plots comparable to figure 12 can be made for other 


responses (e.g., blade material degradation due to oxidation or 
other causes) provided that appropriate factors are used in the 
MFIM. 

The information from figure 12 can be combined with costs 
for fabrication and failure (penalties). The results are shown in 
figure 13 as log-log plots for probability of damage initiation 
versus number of cycles and for total cost versus fatigue cycles. 
Note that the cost increases very rapidly with fatigue cycles 
higher than 10 000. Interestingly, the information in figure 13 
is really the cascading type because it can also be used to 
generate information for benefits accrued by improving mate- 
rial strength or controlling the quality of processing. Costs to 
improve structural reliability by decreasing scatter are more 
beneficial than costs to increase strength for comparable prob- 
abilities in general (unpublished in-house results by the author 
and his collaborators). 
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(b) 


Random variables 


Figure 15. — System reliability using fault tree simulation Figure 15. — Concluded, (b) Random variables. Creep 

and primitive variables sensitivities, (a) Failure mode. coefficient, B 0 . 

(b) Random variables. Creep coefficient, B 0 . 


Fault Tree Simulation for System 
Reliability 

Systems usually fail by combinations of multiple failure 
modes. Multiple failure mode reliability is evaluated by a 
formal combination of the probability of failure of each failure 
mode. Traditionally, the formal method for combining those 
probabilities has been the fault tree simulation. Figure 14 
schematically depicts a fault tree simulation for the SSME 
high-pressure blade. The evaluation includes four failure modes, 
the probability of failure for each being determined by proba- 
bilistic structural analysis. The modes are then combined by 
classical probability methods. The details of this process are 
described in reference 10. Herein, we present some typical 
results and discuss their significance to material behavior and 
its influence on structural system reliability. The parameters 
(primitive variables) with their respective scatter standard 
deviation and assumed distributions included in the simulation 
are listed in table m. Results obtained for the probability of 
system failure are shown in figure 15. The probability for 
system failure from each individual failure mode is shown in 
figure 15(a) with the simulated correlation coefficients for the 
three failure modes. From figure 1 5(a), the system is predicted 
to fail by creep (stress rupture) because that failure mode has the 
highest probability of failure. Also, the stress influences the 
vibration failure mode and the creep failure mode significantly , 
even though the system failure probability from stress only is 
relatively small as compared with the other two. The sensitivity 
factors for system failure probability are shown in figure 1 5(b). 
The dominant factor is the direction solidification angle 6y 


measured from the blade radial (spanwise) axis. The modulus 
and the Poisson’s ratio have about the same influence, but the 
rest have relatively negligible influence. Recommendations for 
materials suppliers and designers are comparable to those 
mentioned previously. Recommendations for the blade manu- 
facturers and the rotor assembler are that they assure that the 
blade solidification direction line up with the blade spanwise 
axis. Important from the previous discussion is that system 
reliability can be formally evaluated for multiple failure modes . 
The critical failure modes are identifiable and their respective 
dominant sensitivities are quantifiable. 

Three approaches were described for evaluating system 
reliabilities. The first is summarized in figure 6, the second in 
figures 11 and 12, and the third in figures 13 and 14, all with 
their attendant discussions. The first and second methods are 
evaluated directly from the probabilistic structural and stress 
analysis, whereas the third requires fault tree evaluation in 
addition to probabilistic structural and stress analyses. The 
computer code described in reference 3 was used for all three 
approaches. 


Concluding Remarks 

A formal method to quantify structural damage tolerance and 
reliability in the presence of a multitude of uncertainties in 
turbine engine components is described. The results of using 
probabilistic methods to predict structural damage tolerance 
and reliability from materials to service environments lead to 
the following remarks: 
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1 . Probabilistic methods via computational simulation can 
be adapted throughout the structural design practice. They 
provide quantifiable information that can be used to reduce the 
costs of product development, certification, and risk. 

2. These methods constitute a virtual statistical desktop 
laboratory applicable at all stages and for all aspects of the 
following: design, material selection and qualification, devel- 
opment, certification, and service life cycles. 

3. Probabilistic methods rely on computational simulation 
results and are primarily useful for decisionmaking, especially 
in the preliminary design stages. 

4. Probabilistic methods provide a significant infrastructure 
that can be used to make predictions for future strategic projec- 
tions and planning to assure better, cheaper, faster products that 
will give competitive advantages. These methods also have 
acceptable reliability, quantifiable risk, and capability to evalu- 
ate the remaining life of existing products. 


Glenn Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, July 25, 1999 
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